Quantum-Mechanical Analysis of Single Particle Level Density* 
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A quantum-mechanical calculation of the single-particle level (s.p.l.) density 
g{e) is carried on by using the connection with the single-particle Green's function. 
The relation between the imaginary part of Green's function and single-particle 
wave functions is used separately for the discrete and continuous states. Within 
the bound-states region the imaginary part of the Green's function is calculated by 
using the wronskian theorem. The Green's function corresponding to the continuum 
is written by using the regular and Jost solutions of the radial Schrodinger equation. 
The smooth part of the rapidly fluctuating s.p.l. density is calculated by means 
of the Strutinski procedure. The continuum component of the s.p.l. density has 
rather close values within either exact quantum-mechanical calculations with the 
Woods-Saxon (WS) potential, or Thomas-Fermi approximation with WS as well as 
finite-square potential wells, provided that the free-gas contribution is subtracted. 
A similar trend is obtained by means of the simple FGM formula for the s.p.l. 
density if the continuum effect is taken into account. 

PACS: 21.10.Ma,21.10.Pc,21.60.Cs 



I. INTRODUCTION 



\ The nuclear level density has a central role in the statistical analysis of nuclear reactions. In- 
formation concerning nuclear densities can be obtained from experiments by means of the analysis 
+j£)f neutron and charged particles resonances, inelastic scattering, and particle evaporation spec- 
tra, so that the theoretical calculation of this quantity is quite useful for validation of nuclear 
^models. Many approaches (e.g. and references therein) have been developed to calculate the 
tingle-particle level (s.p.l.) density on which is based the nuclear level density calculation, its es- 
tablishment being yet considered a difficult task. A constant s.p.l. density g ~ A/14 MeV^ 1 was 
^jised in various statistical model descriptions of the nuclear reactions. An increased criticism has 
^een expressed on this constant value which is inconsistent with the number of A nucleons in the 
nucleus and the usual Fermi energy F = 38MeV. 

The semiclassical approximations could be the starting point of a basic analysis of the s.p.l. 
density, including the energy dependence. First, the Thomas-Fermi approximation proves ade- 
quate for a special class of potential wells [U which describe the nuclear mean field (Woods-Saxon, 
harmonic oscillator, trapezoidal potential well). However, for potentials like square- well and infi- 
nite square potential, the results deviate significantly from the quantum-mechanical calculation. 
Enhanced approximations add the lowest terms in h as corrections to the Thomas-Fermi results 
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An exact quantum-mechanical calculation was provided by Shlomo [[J by using a Green's func- 
tion approach. A particular case is the s.p.l. density for a finite potential well. Since it is including 
the free-gas states, one can calculate and subtract this component by using the Green's function 
associated with the respective single-particle Hamiltonian. Shlomo found by means of both semi- 
classical and quantum methods that, for a realistic finite depth potential well, the s.p.l. density 
decreases with energy in the continuum region (the continuum effect). 

The paper is organised as follows. Section II is an introduction of the general theoretical 
aspects concerning the s.p.l. density, namely the definition, its relation with the Green's function, 
the Strutinski smoothing procedure, and semiclassical methods usually involved. Section III gives 
a general description of the single particle Green's function theory and a method for the calculation 
of its imaginary part. The model parameters for the mean field potential and numerical results are 
discussed in section IV. Finally, the conclusions are drawn in section V. 



II. BASIC FORMALISM FOR SINGLE-PARTICLE LEVEL DENSITY 



A. The quantum- mechanical single-particle level density 



A brief presentation of the s.p.l. density definition, Green's function and relation between them 
are given in this section. The s.p.l. density is defined as [l[ 

g(e) = Tr(8(e - H)) , (1) 

where H is the corresponding single-particle Hamiltonian (mean field) 

H=^ + V(r). (2) 

The eventual bound states of the Hamiltonian given by Eq. (§) can be obtained from the time 
independent Schrodinger equation 

H\na) = e n \na) , 

a being the degeneracy given by other quantum numbers. It is supposed that the normalisation 
for the eigenvectors corresponding to bound states is 

(n'a'\na) = <W<W, (3) 

so that the s.p.l. density following Eq. ([!]) is 

g B {e) = J2 6 ( £ - £ n) ■ (4) 

na 

In the case of a finite potential well which has bound as well as continuous states the s.p.l. 
density can be considered for the two regions of the spectrum 

g{£) = 9b(s) + gc(e) , 

where gsie) and gc(s) are the contributions given by the bound and continuus states respectively. 
The continuum contribution are determined from the scattering phase shifts <%(e) by [III 



9c(e) = -E( 2 i + !)" 
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It can be determined also from Eq. ([[[) but the eigenvectors corresponding to the continuum states 
can not be normalised similarly to the bound states case. The relation @ becoming 

(e'a'\ea) = 5(e — e')5(a — a') , (5) 

it results 

9c(e) = £ / d ^( £ " e')(e'a'\eW) , (6) 

a' J 

where the sum after a' can be an integral if this parameter is continuous. For the moment we will 
assume that the states are not degenerate so that the s.p.l. density in the whole energy range has 
the form 

g(e) = £ 5(e - e n ) + [ de'S(e - e'){e'\e') . (7) 

An important quantity is the number of states N(e) with energy bellow e, which will be used 
latter for illustration of the continuum effect, 



re 



N(e) = f de'g(e') . (8) 

J — oo 

By means of the condition 

deg(e) , 

-co 

Eq. (§) can be also used for determination of the Fermi energy sp which is a significant quantity 
for the description of the nucleus ground state or nuclear reactions analysis. 

B. The Green's function approach 

The single-particle Green's function G(r, r'; e) corresponding to the Hamiltonian (§) is defined 
by the equation 

(e-H)G(r,r';e) = 5(r - r') . (9) 

A first assumption has considered the nucleus in a large volume Q (i.e. within an infinite 
potential well), the continuum being discretized. In this case, the spectral representation of the 
Green's function 

G{iy . E)=j: mm (10) 

can be used to obtain 

lim / dr lim (G(r, r'; e + ia) — G(r, r'; e — ia)) 

Q-+0+ 2ni J r'-+r 

= J_\ im y( 1 1 ) = -ys( £ -e n ) } (11) 

2-ki <x->o+^ \e- e n + ia e-e n -ia) V 

where it has been used the well-known relation || 
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lim 



1 



=1= i7i5(x — a) . 



a->o+ x — a ±ia x — a 

By comparison of the Eqs. (Q) and (|i"l"D it results the following relation between the s.p.l. density 
and Green's function (see also B) 



9\?) 



-— lim f dr lim Irn(Gr(r, r'; e + ia)) . 

7T a^0+ J r— >r' 



(12) 



The demonstration can be generalised at this point for the both sides of the spectrum. Without 
consideration of the nucleus in a finite volume Eq. fliCf ) becomes 



, (r\n)(n\r') f (r\e') (e'\r') 



e — E r 



e — s' 



(13) 



while an integral corresponding to the continuum energy range is added similarly to the sum in 
Eq. (|TTD . Next, within this additional term it can be used the relation 



J dr(r|e')(e'|r) = (e'\e') 



so that the respective continuum contribution becomes just the second term in Eq. (|7|). Therefore, 
Eq. (|12"D is proved worthy not only for the bound states but in the general case. 
Next, the dependence of the s.p.l. density by e and r can be expressed as 



g(e, r) = lim ImG(r, r'; e + ia) 

71 Q^0+ 



(14) 



It is reduced to the dependence by e and r in the common case of a central potential. 

The s.p.l. density for the free particles without spin, i.e. for V(r) = in Eq. (g) can be also 
obtained for Eqs. (0) and (|i~4"D . The single particle Green's function in this case is 



Gi +) (ry;£] 



lim Go( r ) r ' ; e + i°) 



2m exp(i/;;|R|) 
W 47r|R| 



(15) 



where R = r — r' and k = y 2me / h 2 . 
Introducing Eq. ( p!5|) in Eq. flliD 



qnie, r) = lim Im 

V ' 7T R->o 

2m 

= — — -77 lim 



2m exp(ikR) 
_~ h 2 AttR 
sin(kR) 
R 



(16) 



it results the expected characteristic of go(e, r) to be a function of only e, for the case of free 
particles 
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1 /2m\ 3 / 2 



(17) 



Furthermore, by using Eq. (|T^) in Eq. (12) one obtains the well-known formula of the s.p.l. 
density for free particles without spin restricted to a finite volume Q M 
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(18) 



We have to pay attention at the numerical calculation of the limit a — > 0+ in Eq. (|12|) . The 
s.p.l. density could be calculated by smearing of the Green's function with a Lorentz function of 
small width 7. By using the relation 

'V r°o Ac' 

lir ? 1 / t, — vr— r '; £> + io = G ( r ' r '; £ + *7) > (19) 

a-»0+ 7T J-oo (£' — + 7^ 

the smeared level density <7 7 (e) can be obtained from Eq. (|l^) while a is replaced by a finite 
but small value of 7. This parameter should have a value much smaller than hut, the spacing 
between major shells. The method is used in several studies [|]|||7| and has the advantage that the 
knowledge of the spectrum of H is not any more necessary. 

Alternatively, an analytical calculation of the single-particle Green's function is adopted in the 
present work, in order to obtain an exact form of the s.p.l. density within the continuum region 
too. Because of the differences between the regular and irregular solutions of the radial Schrodinger 
equation for the two regions of the spectrum, the respective limits a — > + in Eq. (|TJj) have to be 
calculated separately as it is shown in section III. 



C. The smoothing procedure 

The level density g[e) varies rapidly with energy, due to the shell effects. However, besides its 
rapidly fluctuating part, g[e) contains a smooth part describing the average behaviour of the level 
distribution |||| 

g(e) = g s (e) + Sg(e) . 

For the discrete states, the smooth part of the level density can be obtained from Eq. (O) by 
smearing out each 5 function to a Gaussian of width T H 



2 

-n 



S(e - e n ) -> — ^exp 

where T should be at least of the order of the shell spacing. Although such a Gaussian preserves the 
normalisation of a 5-function, curvature corrections have to be introduced in such a way that not 
only g s (e) but also its first 2M derivatives are reproduced correctly @||. The unitary treatment of 
the bound and continuous states is realized by means of an average of the s.p.l. density (Strutinski's 
procedure) [Q 

/oo 
de'F(e - e')g{e') , (20) 
-00 

where 

Fix) = -^e-(^) 2 4(x 2 /r 2 ) , 

1 y/TT 

1 

and L| f is the associated Laguerre polynomial [[|. The smoothing parameter F has to be chosen 
MM greater than Tluj. 



5 



D. Subtracting the free gas contribution 

The background level density due to the free Fermi gas (which goes like \Je) exists whether or 
not the potential well is finite Q, the former case yielding the quantum fluctuation of the density. 
A subtracting procedure was introduced by Shlomo |l]-f|, the states corresponding to the free 
Hamiltonian H = P /2m being taken away for e > 

g{e) = J dr(0(e,r)-0 O (e,r)) . (21) 

Therefore, the s.p.l. density decreases with energy in the continuum region (the so called continuum 
effect). 

A particular problem occurs for protons. The corresponding free gas contribution has to be 
considered with the Coulomb interaction included. The smoothing procedure (|20"D will be involved 
following the subtraction of the free gas contribution. 



E. Semiclassical calculation of g s (e) 

The semiclassical approximations for calculation of the s.p.l. density are commonly used 
HI [3,10) ■ The ^-expansion (Wigner-Kirkwood) of the single-particle partition function is used 



in this respect [jHj], performing the so-called Wigner transform. The terms with odd power of h 
vanish for the smooth single-particle potentials. Therefore, the Thomas-Fermi (TF) approximation 
(h° term) and the semiclassical (SC) approximation (TF+h 2 correction) are generally considered. 
The numerical calculations in this work have been carried out for the TF-method only. The results 
of the TF approximation by taking into account the spin degeneracy and neglecting the spin-orbit 
interactions are given in several papers |Tl,p|JlO[ 



9tf(s) = ± /dr (e - V(v)f 2 Q(e - V(r)) . (22) 

For the finite well potentials we have to subtract Jl],|2j the contribution of free Fermi gas when 
e > 0, which is given by Eq. (|17D. 



III. SINGLE-PARTICLE GREEN'S FUNCTION CALCULATION 

In this section a partial wave analysis for Green's function is presented for particles with spin 
H. A brief review of some elements which are necessary for such description of particles is given 
firstly. 

Let xt an d C(lsj,muM) be the normalised eigenvectors of the spin, and respectively the 
Clebsch-Gordan coefficients. By using spherical harmonics ^ m (?) there are defined the functions 

i 

Yj£(?) = E C{lsj, muM)Yr(i)xt ■ (23) 

mv 

They are eigenf unctions of the total angular-momentum operator, its z component (with respect 
to a given z axis), the orbital angle- momentum operator, and the spin operator. 
Following the expansion of the Green's function |J 
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9m 1 

G^(r,r>;e) = ™ — £ Y^Y^r')^ (r, r'; e) , (24) 
the partial-wave Green's function G\f\r,r'; e) satisfies the equation 



T i.i 

d 2 1(1 + 1) , 2 2m 



-r —T- ~k 2 + ^V l3 (r) j Gf\r,r';e) = -5(r - r>) (25) 

with the boundary conditions specifying that it is zero at r = 0, for fixed r' and at r — ► oo 
it contains outgoing (for G^) or incoming waves (for G(->). The solution of Eq. @ can be 
expressed by using the regular uij(r) and irregular solutions v\f\r) of the radial equation, 



J ij 

r « v "' ,c;_ w! 



Gh >(r,r;e) = —f , (26) 



'ij 

where W\j is the wronskian of the two solutions, and r < and r> are the lesser and the greater of r 
and r', respectively. The regular solutions are defined for r — > so that it is satisfied the boundary 
condition 

uij(r = 0) = . (27) 

Because of the hermiticity of the Hamiltonian, the regular solution can be chosen real. The irregular 
solutions are defined for r — > oo by 



if e < 

e ±ikr if e > , 



where k = J—2me/h . 



A. Bound-states region (e < 0) 



The imaginary part of G(r, r; e + iot) is obtained following the calculation of the limit a — > 0+. 
Moreover, we have to pay attention to the both possibilities for e to be or not an eigenvalue. In the 
first case, since the regular and irregular solutions can be chosen real and are linearly independent, 
the limit in Eq. (O) is easy to perform and the imaginary part of Green's function becomes zero. 

In the second case the regular u(r) and irregular solution v(r) are linearly dependent, so that 
it results 



uiAr) = Cy%(r) , (29) 



where the coefficient C\j can be obtained from the boundary conditions. 
Let Eq be an eigenvalue, and K a complex number given by 



2m(e + ia) . , 
K o = \l ^2 = m - a , 



where a' — > 0+. The wronskian can be expanded around a' 
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W (K 



dW 



a' + 0(a 



r2\ 



(30) 



K =IK 



while the use of the wronskian theorem leads to J5.ll 



dW 



~ / druj(r) . 
Go Jo 



(31) 



The imaginary part of the Green's function can be calculated now by using Eqs. (|26"D and (j31|). 
The result is a sum of 5-functions 



Im (G(r, r; e)) = D n u 2 n (r)5(e - e n ) . 



(32) 



Finally, by introducing Eq. ( |3"2"D in Eq. (|T2) and comparing with Eq. (ffl), the coefficients D n can 
be obtained and the s.p.l. density at negative energies becomes 

9(£,r) = ^2ul(r)5(£ - e n ) , (33) 

n 

where it has been assumed the following normalisation condition for the regular solutions 

drw^(r) = 1 . 

The solutions of the radial Schrodinger equation considered in this frame 



h d 2 u n (r) 
2m dr 2 



+ 



V(r) 



h 2 1(1 + 1) 

2m r 2 



u n (r) = 



(34) 



have been obtained for V(r) being the sum of the Woods-Saxon (WS) nuclear potential and the 
Coulomb potential in the case of protons. The former has the expression 



V ws (r) 



To 



1 + exp[(r - R)/a ] 



(35) 



where Vq < 0, R and a Q are the depth, the radius and the diffuseness parameters of the well, 
respectively (given within the following section), while the latter is 



VJr) 



2R C ( 3 B? c ) r < Rc 

Z,J r>R r 



(36) 



for a uniformly charged sphere of radius R c . The boundary conditions satisfied by these solutions 
are 



U n (Q) = 



«n( r )lr-oo • 

The eigenvalues of the WS potential well, which can not be found analytically, have been 
obtained by using the Nilsson-orbits computer code developed by Hird fT2fl . A basis of harmonic 
oscillator eigenstates 
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RnAr) 



2A 3 / 2 n! 



r(n + i 



(Ar 2 )^(Ar 2 )exp(-^ 



has been used, with the condition for the best convergence of the WS states in the harmonic 
potential given by the relation [T2I 



A = 0.15466 



Vol 



R 



(37) 



between the size parameter and the parameters of the well, where the energies are in MeV and the 
lengths in fm. Illustrative results of the numerical integration of Eq. fl34|) by using this method 
are shown in Fig. 1. 
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FIG. 1. (a) The Woods-Saxon nuclear potential shape for the global parameter set given, and (b) 
corresponding radial wave functions for neutrons. 



B. Continuum-energy region (e > 0) 

The methods used in this work for evaluation of the s.p.l. density within the bound-states and 
continuum regions, respectively, are distinct. There are two linearly independent solutions for the 
radial Schrodinger equation, defined by various boundary conditions 

Uij(r) — > sin (kr — y + 5^ when r — > oo, (38) 

and 

uij{r) — > cos (kr — l -f + %J when r — > oo . (39) 

Eq. (|38f) is providing the regular solution f5],p|l, while Eqs. (|38|) and ( ]39"D can be used to obtain 
the irregular solution (Jost solution) uL- (r) satisfying the condition given by Eq. (|28"D , 
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v if\ r ) = v>ij(r) ± iuij(r) . (40) 

The wronskian of the two solutions is r-independent (according to the wronskian theorem) with 
the result Wij = —k. Consequently, the imaginary part of Green's function becomes 

bn(^(r,r; e )) = -^. (41) 

Therefore, the r- and e-dependent s.p.l. density for e > is 

9m 1 

9^r) = — 2 -Y,(2j + l)ul(r), (42) 

i.e. it is given only by the regular solution of the radial Schrodinger equation corresponding to the 
potential well. 

The latest equations can be used for the calculation of the s.p.l. density for free particles too. 
The regular solution in this case is given by 

u^(r) = krji(kr) , (43) 

where ji(kr) is the spherical Bessel function, regular in origin, and go(e,r) becomes 

2?7T< 

g (e,r) = ^kr 2 J2^J + l)j?(kr) . (44) 

Since 2j + 1 = 2(2/ + 1) in the general case of a single-particle potential which neglect the spin-orbit 
interaction, it results by using the relation || 

oo 

£(« + !)#(*) = 1, 

1=0 

that go(e,r) is the same as given by Eq. flT7|) except the spin degeneracy is taking into account 

t \ 4m ; 2 

g (e,r) = —^kr . 
Tin 

This result is an useful test for the correctness of the relation fl4"2] ) which has been in thiswork the 
starting point for the numerical calculation of the s.p.l. density in the continuum region. 

The numerical solutions of the Schrodinger equation have been obtained by using the Cowell 
method |13|]14}| . A modified version of the optical model code SCAT2 |14]] has been used in this 
respect. The results of the numerical integration for neutrons is shown in Fig. 2. 

The solutions are conveniently defined by the condition that they are given asymptotically, e.g. 
for neutrons, by 

Uij(r) = sin ^kr — — + when r — > oo , (45) 

where 5ij is the phase shift corresponding to the mean-field potential. In particular, the phase shift 
can be used in the calculation of the differential cross section |5|JTT[1 for the particles scattered on 
a finite-range potential (for the Coulomb potential can be also made). The s-wave phase shift for 
the WS potential well is shown as a function of energy in Fig. 2(b). 
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FIG. 2. The s-wave radial solution of the Schrodinger equation in continuum, for the Woods-Saxon 
potential (solid curve) as well as without it (dotted curve), and (b) the corresponding s-wave neutron 
phase shift, for the nucleus 93 Nb. 



IV. MODEL PARAMETERS AND RESULTS 

The quantum-mechanical and TF calculations of the s.p.l. density for neutrons within the 
nucleus 56 Fe have been carried out by using the following potential wells. The respective parameters 
were determined by looking for the available experimental data on nuclear radii and separation 
energies to be reproduced by the corresponding single-particle Hamiltonian |J. 



A. Potential well parameters 

(i) The harmonic- os dilator (HO) potential as sample of an infinite potential well, but with a 
smooth surface and applicable only for A <100, is given by 

V HO (r) = ~mu; 2 r 2 , (46) 



where the parameterisation 



has been adopted. On the other hand, the usual size parameter of the potential A = muj/h has 
been replaced according to Eq. (|37|) for the basis of HO eigenstates which is used for the expansion 
of the WS potential eigenstates, for a better convergence of the numerical method [12 . 
(ii) The infinite square-well (SQ) potential of the form 

has been used with the depth and radius parameters 
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N — Z 

V Q = -46 + 33— —t 3 (MeV), 

R = R v = 1.12A 1/3 + 1.6 fm, 

Actually, it is the basis for the Fermi-gas model (FGM) of the nucleus. The respective eigenvalues 
are provided by the boundary condition for the single-particle wave functions to vanish at the 
nuclear surface (r = R) while the square-root energy dependence of the s.p.l. density is given by, 
e.g., the semiclassical Thomas- Fermi formula (|22|) . 
(Hi) The finite square (FSQ) potential of the form 



has been used, with the depth 



N — Z 

V = -47 + 33 — —t 3 (MeV) 



where £3=1 for a neutron and -1 for a proton, while the radius has the same value as for the WS 
potential well. 

(iv) The WS potential well parameters within Eq. fl35|) as given by Shlomo's global set are |1[] 

N — Z 

V = -54 + 33 — - — £ 3 (MeV) 

R= fTXT^TpWai ' ^ = 1.12^/3 + 1.0 fm (50) 

[1 + (na/R) 1 '*] 

a = 0.70 fm, 

where R is determined by iteration. 



B. Results and discussion 

The number of the bound single-neutron states with the energy less than e, i.e. N(e) given by 
Eq. (||), is shown in Fig. 3(a) as a function of e for the nucleus 40 Ca, by using the WS potential 
within the both QM and TF methods. It is found similarly to Shlomo that, in the case of the 
WS smooth potential, the TF formula may be considered a good approximation of the quantum- 
mechanical results. Additionally, in Fig. 3(b) is shown the effect of the free-gas contribution 
subtracting, in the frame of the TF approximation. 

The continuum effect on the s.p.l. density for neutrons is firstly shown within the TF method, 
in the case of the nucleus 56 Fe. The use of the infinite square-well provides the usual s.p.l. density of 
the Fermi-gas model. At the same time the respective finite well determines the continuum region 
of the spectrum where the consideration of the free-gas contribution by Eq. (|2l|) is illustrated in 
Fig. 4(a). Next, the same effect in the TF approximation but using the realistic WS potential well 
is shown in Fig. 4(b). It could be noted that, in spite of the quite distinct energy dependences 
due to the SQ and WS potentials wells, the continuum components are rather similar following the 
proper consideration of the free-gas contribution. 

The analysis of the continuum effect is completed by the quantum-mechanical calculation for 
the same neutron level density of the nucleus 56 Fe in the Woods-Saxon potential well. The smooth 
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e (MeV) £ ( MeV ) 

FIG. 3. The number of single-neutron states N(e) as a function of e for the nucleus 40 Ca and the 
Woods-Saxon potential, given by (a) the quantum-mechanical calculated single-particle energies (his- 
togram), and the TF approximation (dotted curve), and (b) the latter method with (solid curve) and 
without (dotted curve) the free-gas contribution subtracted. 




£ (MeV) e (MeV) 

FIG. 4. The s.p.l. density for neutrons of the nucleus 56 Fe as a function of the single-particle energy, 
given by the TF approximation with (a) the infinite square- well (dashed curve), finite square potential 
well (solid curve), and (b) the Woods-Saxon potential (dashed curve), as well as by the smoothed quan- 
tum-mechanical calculation with the same WS potential well (solid curve). The dotted curve correspond 
to the QM calculation carried out without subtraction of the free-gas contribution. 

level density obtained by using the smooth parameter V = 1.2hu, where hu is given by Eq. ([17]), 
and the order of Laguerre polynomial M=2, is shown in Fig. 4(b) too. The comparison with the 
results of the TF approximation, which involves the same potential well, proves the suitability 
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FIG. 5. The same as Fig. 4, for the quantum-mechanical calculation with the WS potential (solid 
curve), TF approximation with the FSQ potential well (dashed curve), and the FGM formula with the 
Fermi energy values of -F=38 MeV (dotted curve) and F=20 MeV (dash-dotted curve). 

of the latter similarly to that previously by Shlomo [|TJ. 

Actually, it results that the continuum component of the s.p.l. density has rather close values 
within either exact QM calculations with the WS potential, or TF approximation with WS as 
well as FSQ potential wells, provided that the free-gas contribution is subtracted. Moreover, 
a similar trend is obtained by means of the simple FGM formula for the s.p.l. density (e.g., 
|T5|] ) , if the continuum effect is taken into account (Fig. 5). The level of agreement between this 
phenomenological formula and the QM calculation could be improved for lower values of the Fermi 
energy F which is usually considered 
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It is thus even possible to reproduce the 
quantum-mechanical results by using the average Fermi energy F ~20 MeV which is determined 
by the nuclear-surface localization of the first nucleon-nucleon collision in preequilibrium reactions 
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FIG. 6. The smooth s.p.l. density for the nucleus Ca as a function of the single-particle energy, 
given by (a) quantum-mechanical (solid curve) and TF approximation (dashed curve), for neutrons, and 
(b) QM calculations for neutrons (dotted curve) as well as protons (dashed curve) and their sum (solid 
curve) . 

Finally, the smoothed results of the quantum-mechanical calculation for the s.p.l. density of 
neutrons as well as protons in the nucleus 40 Ca are shown in Fig. 6(b). In the proton case the 
calculation of the free-gas contribution has included the Coulomb interaction too. This case was 
chosen especially due to the equal numbers of neutrons and protons of 40 Ca, so that it can be 
observed the effect of taking into account the both kinds of nucleons. Actually, the comparison of 
the exact QM and TF approximation is shown in Fig. 5(a) only for neutrons, while next is given 
the QM calculated contributions of the both systems of nucleons as well as their sum. 

V. CONCLUSIONS 

A quantum-mechanical calculation of the single-particle level (s.p.l.) density g(e) is carried 
on by using the connection with the single-particle Green's function. The relation between the 
imaginary part of Green's function and single-particle wave functions is used separately for the 
discrete and continuous states. Within the bound-states region the imaginary part of the Green's 
function is calculated by using the wronskian theorem. The Green's function corresponding to the 
continuum is written by using the regular and Jost solutions of the radial Schrodinger equation. 
The smooth part of the rapidly fluctuating s.p.l. density is calculated by means of the Strutinski 
procedure. 

Recently Shlomo || calculated the s.p.l. density by using the Green's function method, the 
phase-shift approach, and the TF approximation. The following different points could be underlined 
between the present work and the Shlomo's analysis. First, the smearing procedure used by Shlomo 
is replaced with the exact calculation of the imaginary part of the Green's function carried 
out distinctively for the both bound and continuum regions of the spectrum. Second, Shlomo 
has considered the nucleus contained in a large but finite box within both the Green's function 
method and the phase-shift approach. This condition is not required in the present work. Third, 
we have taken into account the case of the protons too, where the Coulomb interaction has to be 
also included in the contribution of the free gas. 

Finally, it results that the continuum component of the s.p.l. density has rather close values 
within either exact QM calculations with the WS potential, or TF approximation with WS as well 
as FSQ potential wells, provided that the free-gas contribution is subtracted. A similar trend is 
obtained by means of the simple FGM formula for the s.p.l. density, if the continuum effect is 
taken into account. The level of agreement between this phenomenological formula and the QM 
calculation could be improved for lower values of the Fermi energy F which is usually considered 
~38 MeV [0]. It is thus even possible to reproduce the quantum-mechanical results by using the 
average Fermi energy F ~20 MeV [fT/|-|19fl. Additional analysis is yet required by, e.g., adoption 
of realistic mean-field potentials energy dependent ||, and systematic calculations of s.p.l. density 
including the continuum effect. 
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